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Abstract. Excitation energies to isobaric states, that are analogs of ground states, are dominated by 
contributions from the symmetry energy. This opens up a possibiHty of investigating the symmetry 
energy on nucleus-by-nucleus basis. Upon correcting energies of measured nuclear levels for shell 
and pairing effects, we find that the lowest energies for a given isospin rise in proportion to the 
square of isospin, allowing for an interpretation of the coefficient of proportionality in terms of a 
symmetry coefficient for a given nucleus. In the (A,Z) regions where there are enough data, we 
demonstrate a Z-independence of that coefficient. We further concentrate on the A-dependence of 
the coefficient, in order to learn about the density dependence of symmetry energy in uniform matter, 
given the changes of the density in the surface region. In parallel to the analysis of data, we carry out 
an analysis of the coefficient for nuclei calculated within the Skyrme-Hartree-Fock (SHF) approach, 
with known symmetry energy for uniform matter. While the data from isobaric analog states suggest 
a simple interpretation for the A-dependent symmetry coefficient, in terms of the surface and volume 
symmetry coefficients, the SHF results point to a more complicated situation within the isovector 
sector than in the isoscalar, with much stronger curvature effects in the first. We exploit the SHF 
results in estimating the curvature contributions to the symmetry coefficient. That assessment is 
hampered by instabilities of common Skyrme parameterizations of nuclear interactions. 
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SYMMETRY ENERGY 

Symmetry energy describes how bulk nuclear energy evolves with changes in neutron- 
proton asymmetry. It is naturally of great interest in the context of extrapolations from 
nuclear systems observed under laboratory conditions to neutron matter. Moreover, 
because of different contributions of the two- and few-nucleon interactions to neutron- 
proton symmetric and to neutron-dominated systems, it is of interest in sorting out the 
impact of those interactions on nuclear properties. 

Energy of a nucleus with neutrons and Z protons, N -\-Z = A, may be represented 
in the form 

E{N,Z) = Eo{A) +Ei+Ec + E^, . (1) 

Here, Eq and Ei stand for bulk nuclear contributions to the energy, with Eq - the 
energy of a symmetric system and Ei - the symmetry energy. Further, Ec is Coulomb 
energy and fimic represents microscopic corrections including shell effects and pairing. 
Charge symmetry implies that the symmetry energy should be a quadratic function of 



the asymmetry {N — Z): 

E^ = ^iN-Z)\ (2) 

Here, aa{A) is a generalized mass-dependent (a)symmetry coefficient [1]. Charge sym- 
metry considerations for uniform nuclear matter, with neutron and proton densities of 
Pn and Pp, respectively, p = Pn + Pp, and Coulomb interactions switched off, similarly 
yield: 

f (P»,Pp) = + ^ = §^(P)+S(P) i'^Y- (3) 

In the literature, both Ei and the coefficient ^(p) are referred to as symmetry energy. The 
density-dependent symmetry energy S is usually expanded around normal density po: 

S{p) = aa+:z— — + •••• (4) 
3 Po 

On account of energy of symmetric matter Eq minimizing at po, at nuclear densities, 
pressure in neutron matter with energy E = Ei+S, sustaining neutron stars, tends to be 
dominated by the symmetry energy. Naturally that pressure is strongly sensitive to the 
value of L in (4). We represent the value of the symmetry energy at po with a symmetry 
coefficient in (4), S{po) = a^, because the generalized symmetry coefficient of Eq. (2) 
should reduce to 5(po) in the limit of A — > oo. 

Changes in nuclear density across the surface, interplaying with the density depen- 
dence of the symmetry energy in uniform matter, should generally make the symme- 
try coefficient dependent on mass. In fact, a simple consideration of the effects of 
volume-surface competition on symmetry energy [1,2] produces the result 

+ ^ (5) 
fla(A) al 

As A — > oo, the coefficient aa{A) reduces here to a^; of is the surface symmetry coeffi- 
cient related to the p -dependence of S and to L in particular [1]. Following the above, 
establishing the mass dependence of the coefficient could provide information on S{p) 
and constrain pressures in neutron stars. 

Unfortunately, when trying to deduce an A-dependence of the symmetry coefficient, 
by fitting a bulk-energy formula with (2) to observed nuclear masses, the results turn 
out to be ambiguous [2], even for such a simple parameterization as in Eq. (5). This is 
because the effects of the symmetry energy compete with the effects of Coulomb energy 
and even with terms that contain no asymmetry, due to the mass-asymmetry correlation 
induced by the line of stability. Fortunately, extending the bulk-energy formula, to reach 
consistency with the charge invariance of nuclear interactions, allows for an analysis of 
the symmetry energy on a nucleus-by-nucleus basis, independently of Coulomb energy 
and energy of symmetric matter. 



ANALYSIS OF ISOBARIC ANALOG STATES 



Extension of the bulk-energy formula relies on the observation that the symmetry en- 
ergy (2) may be rewritten in terms of isospin as 

£i = l^r2 (6) 

A 

However, nuclear contributions to the energy should be scalars under rotations in isospin 
space [3], with z and transverse isospin components treated democratically. Thus, the 
symmetry energy should be rather written as 

= — ^^^2^ — e}U.(Tf + Ti)= — ^^T(T+l). (7) 
A A ^ ^ ^ A 

With this extension for the symmetry energy, one can consider states of lowest energy in 
a given nucleus, for a fixed isospin T, rather than absolute ground states that normally 
have T = |r^|. Those excited states end up being isobaric analog states (IAS) of ground 
states of other nuclei within the same isobaric chain. Regarding the bulk contributions 
to the energy, the N — Z portion is the same for those states as for the ground state of 
the specific nucleus and the Coulomb contribution is the same as well. The excitation 
energy to such a state should then reflect change in the symmetry energy associated with 
the change in isospin 

^lAS = ^^AT'+AEr^. = 1^ (ri - |r,|) +AE^, , (8) 

where the isospin in the last term pertains to IAS and where we used T = \Tz\ for the 

ground state. 

To test whether Eq. (8) could be valid, we use the measured energies of IAS, primarily 
from the compilation [4], and construct excitation energies corrected for microscopic 
and deformation effects, following the results of [5]. For selected narrow intervals in 
mass, we plot those energies in Fig. 1, vs scaled difference in isospin squared, between 
the IAS and ground state. If the generalized symmetry coefficient varies slowly with 
mass, the excitation energies should line up, according to Eq. (8), along the lines passing 
through the origin, with a slope given by the symmetry coefficient. It is seen in Fig. 1 
that the excitation energies indeed line up in such a manner, with line slopes increasing 
as the mass number increases. 

The generalized symmetry coefficients, obtained by fitting the IAS excitation spectra 
at a given A, are next shown as a function of A in Fig. 2. The errors in the figure 
are estimated assuming residual unsubtracted microscopic contributions to the input 
energies, of an r.m.s. magnitude of d^ms = 0.5 MeV. A significant variation with mass 
is observed for coefficient values. The typical expectation of > 20 MeV is met for 
mass numbers A > ICQ. However, for A < 10 even as low values as ~ 10 MeV are 
reached. It is apparent that for masses A > 20, the dependence of the coefficient on mass 
is well described by the volume- surface competition formula of Eq. (5). Depreciation of 
coefficient values with mass, at lower mass values, is faster than predicted by (5). 



T 1 1 1 1 1 1 1 1 1 1 1 1 r 




12 3 



4(t!-|TJ)/A 

FIGURE 1. Excitation energies to the ground-state IAS, corrected for microscopic effects [5], plotted 
vs change in T{T + 1) scaled with A, represented by different symbols for different indicated intervals in 
the mass number A. The solid lines represent linear least-square fits to the represented energies. The fits 
are forced to pass through the coordinate origin. The dashed vertical line helps to read off the values of 
aa{A) for the specific A-intervals, from intersection of the fitted lines with the vertical. 
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FIGURE 2. Generalized mass-dependent asymmetry coefficient a„(A) extracted from excitation ener- 
gies to ground-state IAS within individual isobaric chains A, when applying shell corrections by Koura et 
a/. The line represents a fit at A > 22 assuming a combination of volume and surface symmetry terms. 
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FIGURE 3. Generalized symmetry coefficient for individual nuclei, within selected isobar chains char- 
acterized by mass numbers A, presented here as a function of asymmetry. The horizontal dash lines 
represent the coefficient values for the isobaric chains as a whole. 



We further test whether there is evidence for a Z-dependence of the generalized 
symmetry coefficient, in the context of claims made in Refs. [6, 7, 8]. Exemplary values 
of coefficients for different members of isobaric chains, with emphasis on the span of 
masses and coefficient values, as well as the number of members in a chain that could 
be analyzed, are shown in Fig. 3. We find no evidence for the variation of coefficients 
with Z. 



FROM IAS TO INFINITE MATTER 

Given the success of the fit with (5), to the results from IAS displayed in Fig. 2, it would 
be natural to proceed directly to the conclusions on properties of infinite matter [9]. Thus, 
the coefficient ~ 35.0MeV could be identified with the value of the symmetry energy 
at normal density S{po). Moreover, the surface symmetry coefficient lO.OMeV 
could be identified with that for the half- infinite nuclear matter [1]. From correlation 
between af and slope L of symmetry energy, within the calculations of half-infinite 
matter [1], one could then deduce L ~ 95MeV. A more cautious strategy involves 
testing of the procedure within a semi-realistic description of nuclear properties such 
as Skyrme-Hartree-Fock (SHF) [10, 11]. 

The complication with SHF is in the fact that the approach violates charge invari- 
ance, meaning that the procedure such as for the data cannot be applied. Another com- 
plication is in the lack of shell corrections established in such a manner as for the 
data. Correspondingly a different way of arriving at generalized asymmetry coefficients 
needs to be found [11]. When the coefficients are found, it is observed that for virtu- 
ally every Skyrme interaction the approach to the limit of clear surface- volume sepa- 



ration is much slower for the symmetry energy than for the energy of symmetric mat- 
ter. In consequence, for nuclei occurring in nature, the coefficients in Eq. (5) turn out 
to be effective rather than representing nuclear matter. Interestingly, for many Skyrme 
parameterizations in the literature, the limit of Eq. (5) with coefficients truly repre- 
sented nuclear matter is actually never reached, even for nuclei calculated [10] using 
A ~ (100000 — 1 000000), while switching off Coulomb interactions. This is because of 
instabilities [12] for many of the interactions, particularly iso vector in nature, that come 
into play in large systems. 

CONCLUSIONS 

Symmetry coefficients can be obtained for individual nuclei, exploiting charge invari- 
ance of nuclear interactions. The coefficients rise with mass number within the range of 
~ (10 — 23) MeV. For A > 20, the coefficients are well described with a surface- volume 
formula a-'^{A) ~ (al)-^ + {al)-^A-'^/^, where ~ 35.0MeV and ~ lO.OMeV. 
Analysis done with spherical SHE calculations suggests that the volume and surface co- 
efficients from the fit are effective in nature rather than representing the infinite matter. 
Assessment of the properties of the latter, using IAS results, is still in progress. 
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